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Abstract. We present here a self-consistent, tridimen- 
sional model of a disc galaxy composed by a number of 
ellipsoidal distributions of matter having different flatten- 
ing and density profile. The model is self-consistent and 
takes into account the observed luminosity distribution, 
the flattening profile and the stellar rotation- and velocity 
dispersion- curves. In this paper we considered the par- 
ticular case of a disc galaxy composed by two spheroidal 
bodies: an exponential disc and a bulge following the r^^^ 
law. 

We studied the behavior of the stellar rotation- and 
velocity dispersion- profiles along the sequence of SOs and 
Spirals, identified by an increasing disc-to-bulge ratio. In- 
side every class, kinematic curves were produced by chang- 
ing the relative concentration of the two components and 
the inclination of the galaxy with respect to the line of 
sight. The comparison with observational data requires 
only two scaling factors: the total mass of the galaxy, and 
the effective radius. 

The model allows also to detect the presence of 
anisotropy in the velocity distribution. In the special case 
of SOs, we explored the sensitivity of the kinematics of the 
model by changing the anisotropy and the flattening of 
the bulge. For intermediate flattening (0.4 < b/a < 0.85) 
it is possible to distinguish a change of anisotropy of 15% 

To show a real case, the model has been applied to the 
photometric and kinematic data of NGC 5866. We plan 
to apply these models to a larger database of SO galaxies 
in a future paper. 

Key words: Galaxies: fundamental parameters - Galax- 
ies: general - Galaxies: kinematics and dynamics - Galax- 
ies: structure - Galaxies: individual: NGC 5866 
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1. Introduction. 

The study of the light distribution and the stellar kine- 
matics in disk galaxies is an important tool to understand 
their structure. It is the only way to study stellar systems 
which are gas poor such as the SOs, and may be crucial 
to distinguish between different scenarios of galaxy forma- 
tion. It allows to discuss if a bulge is flattened by rotation 
or by anisotropic residual velocities and how the different 
components of a galaxy mutually interact. 

In the past, in front of the success achieved in study- 
ing the gas kinematics of disk galaxies, much less attention 
has been paid to the modeling of their stellar kinematics 
from the observations. A reason of this was the lower ex- 
tension of the available stellar rotation curves, based on 
absorption lines, with respect to that obtained with emis- 
sion lines (optical or 21cm). However, the inner part of 
the galaxies is often poor of cold gas, and can be studied 
with data from the stellar component only. In addition, 
the progresses in detectors and data analysis techniques 
(Sargent et al. 1977, Bertola et al. 1984, Kuijken & Mer- 
rifield 1993) extended the range explored with stellar data 
and produced a large quantity of stellar rotation- and ve- 
locity dispersion- curves, making compulsory the creation 
of realistic theoretical models for their interpretation. 

However, the modeling of stellar components in disk 
galaxies must overcome several difficulties. First, the non- 
negligible thickness of stellar disks, compared to the gas 
ones, that requires a tridimensional approach. Second, the 
fact that the observed rotation curve does not represent 
the circular rotation defined by the potential, because of 
the partial transparency of the stellar body and of the 
presence of the velocity dispersion. The stellar rotation 
curves are integrated along the line-of-sight across regions 
of different kinematics and are influenced by the gradient 
of mass of these regions. As consequence, the line-of-sight 
velocity distribution present deviation from the classical, 
Gaussian shape. The influence of the velocity dispersion 
and its eventual anisotropy are particularly important in 
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early-type galaxies, where the velocity dispersion is of the 
same order of magnitude of the rotational velocity. This 
last effect is known with the name of asymmetric drift 
and the first attempt to evaluate it from observations date 
back to 1961, with the van der Hulst's work on the obser- 
vations of NGC 4111 made by Humason and Oort (Van 
der Hulst 1961). More recently, the influence of the asym- 
metric drift on the stellar rotation curves has been taken 
into account by some authors ( Bertola fc Capaccioli 1975t 
Kormendy 1984, lUingworth & Schechter 1982, Fillmore 
et al. 1986, Zeilinger et al. 1990). These theoretical mod- 
els and interpretations of stellar kinematics were assuming 
an isotropic velocity dispersion or were limited to edge-on 
galaxies. Most part of models were simulating the galaxy 
by means of the disk component only, with the only excep- 
tions of the Fillmore et al. (1986) model and the Zeilinger 
et al. (1990) work. Both papers were based on non self- 
consistent models and the first one made use of photomet- 
ric data only. 

In the last few years, a number of self-consistent, multi- 
components dynamical models for early-type disk galaxies 
have been presented (e.g. Cretton & van den Bosch l"999| , 
Emsellem et al. 1999; see also Merritt 1999 for a review). 
Some of them are based on the Multi-Gaussian Expan- 
sion approach (Emsellem et al. 1994, Loyer et al. 1998) 
for inferring the mass distribution from the observed 
luminosity; most of them adopt the orbit-superposition 
method ( ^chwarzschild 1979] Emsellem et al. 1994, Rix 
et al. 1997, van der Marel et al. 1998, Cretton & van den 
Bosch 1999, Emsellem et al. 1999) in order to derive a ve- 
locity distribution able to reproduce the observed kine- 
matics (including the higher-order moments). These works 
have established the importance of taking into account 
the effects of the asymmetric drift, the projection along 
the line-of-sight and the deviation from the pure Gaus- 
sian shape of the line-of-sight velocity distribution due to 
the superposition of the different components with differ- 
ent kinematic behaviors. 

While this last approach remains the most general and 
flexible way of deducing the dynamical properties of a sin- 
gle galaxy from the observed photometry and stellar kine- 
matics, its application to a large database of galaxies - in 
order to derive general relations between physical param- 
eters - is not straightforward, due to the large number of 
parameters involved. A model with few parameters to be 
constrained by the observations could be a better choice 
when dealing with a large number of objects. 

We present here a self-consistent, multi-component 
model of disk galaxies based on an analytic approach to 
the problem. We also try to derive the relations between 
the observed properties and the model parameters that 
can be useful to the observers in interpreting the data of 
a single galaxy. The real case of an SO galaxy (NGC 5866) 
is also presented and modeled. 



2. The hypothesis of the model 

2.1. The density profiles 

We want to describe a stellar system similar to a real 
one, and in particular to a disk galaxy. A first step is to 
define a suitable distribution function f{x,v) containing 
the tridimensional structure of the galaxy and its velocity 
field. 

The structure of the disk galaxies can be modeled by 
the mixing of different components: a bulge, a disk and 
eventually other components (lens, nuclear disk, halo). 
Every single component may be reproduced by an ellip- 
soidal density distribution whose projection on the sky 
generates elliptical shapes. If the axial ratio of each one 
of these components is constant with radius, but differ- 
ent for every component, we can apply a generalization of 
the Newton's theorems for spherical bodies, that makes 
simpler the calculation of the kinematic properties. This 
approximation allows anyhow to generate projected light 
distributions with variable ellipticity, as observed. This is 
due to the predominance of the different components at 
different distances from the nucleus. An example of a com- 
posite light distribution - which is rounder in the center 
and in the outer regions, while flattens at intermediate 
radii - is given in Fig. 0. Spindle isophotes may also be 
reproduced in this way. 

Another choice concerns the shape of the different 
components. Differently from the studies on the ellipti- 
cal galaxies, the works made on the shape of the disks 
reveals slightly or absent deviations from the oblateness 
(Sandage et al. 1970, Grosbol 1985, Binney & de Vau- 
couleurs 1981, Magrelli et al. 1992). The mean deviation 
suggested should be < 5%. A different situation is found 
in some bulges, whose structure may range from slight 
triaxial (Bertola et al. 1991) to strongly triaxial in rare 
cases (Stark 1977). 

Taking into account all these facts, we decided to as- 
sume a cylindric symmetry for the galaxy's components. 
This seems to be a good compromise between a "realis- 
tic" description of a single galaxy, yet with a large number 
of free parameters, and a general deduction of the galaxy 
dynamics and structure with the minimum number of con- 
straint obtained from the observations. 

Given the previous hypotheses, the density profile of 
each component can be described in terms of the radial 
coordinate only, with the values of their parameters con- 
strained by observations. 

As shown by Young (1976), the r^/'* luminosity pro- 
file followed by the elliptical galaxies and the bulges can 
be deprojected obtaining an intrinsic density law. In the 
last years, evidence has been presented that the Young law 
can indeed be a poor fit, and the Sersic law, r^/", was pre- 
sented as a better-suited function (Caon et al. 1993). By 
analogy, the case of spiral-galaxy bulges has been revisited 
and new trends have emerged, showing a large spread in 
n peaking at n ~ 1 (exponential bulges), but ranging up 
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Fig. 1. Left: Isophotes of a model composed by a spherical bulge following the R^^* law, and an exponential disc with 
b/a = 0.1. Right: The global axial ratio of the model shown in ^ as function of the radius. Note the peak of isophotes flattening 
at intermediate radii, obtained with two components of constant flattening. 



to n ~ 15). (Andredakis & Sanders 1994; de Jong 1997; 
Andredakis et al. 1995; Broeils & Courteau 1997; Prugnel 
& Simien 1997). The deprojection of the general r^/" law 
can then approximated by the analytical formula (Ciotti 
1991): 

« ^(2n-l)/2» (1) 

where x=r/i?oB and b « 1.99865n — 0.326968, as given by 
Caon et al. (1993), which is valid for every value of the 
flattening of the galaxy. i?cB is the bulge effective radius, 
corresponding to the distance from center in which half of 
the total bulge luminosity is contained. 

For the disk component, the observed exponential pro- 
file can be deprojected generating the law: 

p{r) = poBaKo{ar) (2) 

where Kq is the incomplete Bessel function of order zero, 
and 1/a is the disk scale length. The disk effective radius 
is i?eD=l-6784/Q;. It is important here to note that the 
"disk" component we adopted is not a real, infinitely flat 
disk, but a spheroid of finite thickness and exponential 
profile. 

The above laws may be cut at a radius i?max for com- 
putation purposes, and put together to produce a very 
wide range of complex global density distributions. The 
introduction of an (usually dark) halo among the galaxy 
components is needed when kinematics at large galacto- 
centric radii are considered. In this paper, we shall explore 
structures composed by bright components only, omitting 
the introduction of a dark halo. This choice follows the 



fact that in the most part of stellar rotation curves pre- 
sented in the literature the regions explored are limited to 
the luminous component inside 1 ~ 2i?e. 

2.2. The velocity distribution 

When the spatial dependence of the distribution function 
/(x, v) is determined by the combination of the above den- 
sity laws, the velocity distribution can be obtained by in- 
tegrating the coUisionless Boltzmann equation producing 
the moments equations. This will make the galaxy model 
self-consistent, warranting its dynamical equilibrium, such 
as in real systems. 

In order to close the system of the moment equa- 
tions, some additional assumptions are needed on the 
model, that may be obtained by imposing some restric- 
tions on the shape of the velocity ellipsoid. The most 
widely used assumption in literature at this regard is to 
assume aim — <^zz everywhere in the galaxy, and then 
to choose some easy-to-handle functional form of the de- 
pendence of the azimuthal velocity dispersion u^^ from 
the radius and the radial velocity dispersion (see, e.g., 
Satoh 1980). However, this hypothesis implicitly describes 
a system with two integrals of motion or /(i?, Lz) (Binney 
et al. 1990) that has been shown (Merritt 1999; Emsellem 
et al. 1999; Cretton et al. 1999) to fail in correctly repro- 
ducing the observed behavior of many objects. 

In order to explore the three-integrals distribution 
functions /(i?, L^, I3) one is forced to relax the aim = <^zz 
assumption (see, e.g., Dehnen & Gehrard 1993). The sim- 
plest way to do that is to assume aim — o'4><t> everywhere 
in the galaxy and that the ratio between the vertical and 
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the radial velocity dispersions is constant in the galaxy 
(thus including the special case of isotropic models). 

We point out that with this choice we are presenting 
just a special set of anisotropic models involving a third 
integral. Nevertheless, it is likely that the models can be 
used to derive a "mean" vertical anisotropy for every ob- 
served galaxy, useful to derive general relations in large 
set of galaxies. For a deeper insight of the kinematic struc- 
ture of a single anisotropic galaxy, a more complex model 
is needed. 

More in detail, for every different component, we as- 
sume the following constraints: 



p— az 



ji d\npalj^{R,z) 



'RR' 



dlnR 



^dR 



(5) 
(6) 



The partial derivatives of the potential in these equa- 
tions can be calculated by making use of the relations for 
a single elliptical component and of the additive property 
of the potential: 



The velocity distribution is locally Gaussian. This do 
not imply a Gaussian total velocity distribution; 

The velocity ellipsoid's principal axes remain aligned - — = 2'kGz\J \ — e: 
with the cylindric coordinate system (i?, </>, z)^ so that 
(^Rz = ctr^ = 0-02 ~ 0; this assumption describes a sit- 
uation similar to that of the solar neighborhood, where 
the deviation of the velocity ellipsoids from the cylin- 
dric coordinate axes is low. 

The ratios of the components of the tensor af^ are con- 
stant inside the galaxy, so that we can write 
2 



(r + 1) + 1 - e,' 



2 / Piim^) 



(T + l)(T + l-e2)3A 



dr(7) 



M8) 



We defined, following the notation of Binney & 
Tremaine (1987), 



1 + T r+l-e^ 



(9) 



a- = s{x)a{0) (3) 

where s{x) is a scalar function of the coordinates. The 
choice produces velocity ellipsoids with same shape in 
every point of the galaxy. 
— A final requirement is the adoption of: cr^^ — grr, in 
every region of the galaxy. 

Under these limits, needed to simplify the equations 
and to reduce the number of free parameters, we can de- 
scribe the velocity dispersion field by means of Ufm^R, z) 
and the homogeneous anisotropy parameter (3, where 



f3- 



'RR 



(4) 



Note that we do allow different components of the 
galaxy to have different /? values. For instance, we can 
assume the disk to have a velocity distribution more 
isotropic of the bulge, as well as the contrary, and cal- 
culate the global line-of-sight velocity distribution. Such 
differences are not only expected, but also at least hinted 
by the applications of models using only one single value of 



anisotropy in a galaxy (Cretton & van den Bosch 1999) 



3. The equations 

A second step is to derive the main galaxy parameters to 
be directly compared with the observed photometric and 
kinematic data. 

Independently from the chosen density distribution, 
the Jeans equations in cylindrical coordinates for a steady- 
state, axisymmetric system having a velocity ellipsoid de- 
scribed by the above equations produce: 



where e — ^\ — b/a is the eccentricity of the component. 

In the regions where the bulge and disk luminosities 
are comparable, the superposition of the rapid rotation 
of the disk to the slower rotation of the other compo- 
nents produces in general a non-Gaussian, 2-peak line-of- 
sight velocity distribution (LOSVD). To obtain the pro- 
jected {/, y, cr} for this distribution to be compared with 
the data is not an easy task, since we must integrate the 
whole LOSVD along the line-of-sight. However, the ob- 
served line-of-sight velocity distribution is often parame- 
terized in terms of {V, ct, /13, /14} of a Gauss-Hermite se- 
ries( as described by van der Marel & Franx, 1993; see 
Appendix A). We chose to evaluate the value of these pa- 
rameters, instead of deriving the full LOSVD from the 
unprojected velocity distribution. 

First, we have to correct the V, a unprojected values 
by properly taking into account the effect of anisotropy 
and inclination of the galaxy. The rotational velocity along 
the line of sight is < F >= V smicoscj), where 4> is the 
azimuthal coordinate and i is the angle of inclination of 
the galaxy (z = 90° corresponds to an edge-on galaxy); 
the component of the velocity dispersion in the direction 
of the line of sight can be expressed as cr^ (l — (icos^i) 
so that the second order projected moment is < >= 
V'^ SIT? i cos^ (f> + a'^ (1 — l3 cos^ i) . 

We now integrate along the line of sight the 5 projected 
momenta {p, <V >,<V'^ >,<V^ >,< >}. This can 
be done easily as long as we assume that every component 
has a Gaussian local velocity distribution, by making use 
of the general formula: 



<Q> (i?pi.oj,^) 



p{R,z)Q{R,z)RdR 



2 / Rproj 



(10) 



p{R,z)RdR 



E. Pignatelli & G. Galletta: Multi-component models for disk galaxies. I 



5 



where the notation < Q > (iJproj , z) is used to indicate 
the generic quantity Q integrated along the line of sight at 
the position (i?proj , z) on the plane of the sky. The values 
of {V, cr, /i3,/i4} cannot be obtained by a straightforward 
integration along the line of sight, because of their intrin- 
sic non-linearity; instead, we can integrate the velocity 
moments, which are linear. Also note that we do not re- 
strict our investigation to the case of major-or minor axis 
profiles. 

Then we make use of the projected values of such mo- 
menta to calculate the skewness ^3 and of the kurtosis ^4 of 
the LOSVD at any given radius. The {V, a, /13, /14} values 
can now be obtained with the help of the approximated 
formulas (van der Marel & Eranx 1993): 

-h3 = 6/(4\/3) 

- /i4 = (C4 - 3)/(8V6) 

-a = {<V^>p-<V >2)/(l + /l4^/6) 

-V=<V>p-a-h3-\/3 
4. The pareimeters 

The model described in the previous sections is then de- 
termined by 4n + 1 free parameters, where n is the num- 
ber of components of the model: the total mass M, the 
scale length R^, the eccentricity e and the parameter /3 
defined above, plus the inclination angle of the galaxy. In 
the simple case of a bulge-|-disk model, this lead to 9 free 
parameters. 

Two of these parameters can be considered scale fac- 
tors respectively for the mass (and so, the kinematic quan- 
tities) and the lengths. We want to choose these scale fac- 
tors in a way that makes simpler the comparison with 
the observations, so we select the global mass Mtot of the 
model and the global Re (defined as the major axis isopho- 
tal radius inside which lies half of the total luminosity) as 
the scale factors. 

The models are then presented in a dimensionless form, 
having, in the bulge-|-disk case, seven (dimensionless) pa- 
rameters: the two axial ratios, the ratios ij = Mb/Mu 
and ^ = Res/ReD, the two values of 0b and /3d and the 
ratio of the (M/L) values. The parameter 77 represents 
the 'physical' bulge-to-disk ratio and the ^ the concentra- 
tion of the bulge with respect to the disk. In the general, 
n— component model, we would face 5n — 3 parameters. 

Not every region of the parameters space is interesting; 
on the one hand, as will be shown in next section, for some 
combination of the parameters wide regions of the galaxy 
model cannot reach dynamical equilibrium. On the other 
hand, some of these parameters are never observed in real 
galaxies: in particular, we kept fixed a {b/a)-^ = 0.2 value, 
which seems a reasonable value for the observed disks, and 
we restricted ourselves to the case where the M/ L ratio of 
the bulge component is greater of the corresponding disk 
value. 

In the simplest case of a bulge-|-disk model, we do ex- 
pect that the four remaining free parameters will cover 



two different physical characteristics: the couple (r], ^) will 
somehow reproduce the morphological behavior of galax- 
ies along the Hubble diagram, while the couple (eB,/3) 
will give us informations about the dynamic of the bulge 
( rotation- or pressure- supported ) of the galaxy. 

5. The results 

The model offers several possible applications of the stel- 
lar component of disc galaxies to the observed data. We 
choose here to illustrate two general applications: the de- 
duction of the mass ratio between the disk and the bulge 
and the possibility to detect the presence of anisotropics in 
the velocity dispersion field. In the end, we will discuss in 
detail the case of a well-known SO galaxy, NGC 5866, for 
which both photometric and kinematic data are available. 

5.1. The models and the Hubble diagram 

The classification of galaxies in a single unificating scheme 
has been one of the earlier results in the study of stellar 
systems. The details of this classification are, however, 
still under discussion. Along the sequence of disk galaxies, 
the classification of galaxies may be expressed in terms 
of photometric properties, such as the relative weights of 
the bulge and disk component and their relative concen- 
tration. Our model allow a dynainical deduction of these 
properties. We can compare observed velocity- and veloc- 
ity dispersion- curves with that of a model galaxy, pro- 
jected at the same viewing angle of the real galaxy and 
integrated along the line-of-sight. In such a way, one can 
deduce the disk-to-bulge ratio of the galaxy in terms of 
mass ratio and relative concentration of the two compo- 
nents. We expect that this dynamical description may give 
different results with respect to the disk/bulge ratio pho- 
tometrically deduced, and allow to correct the M/L ratio 
assumed and eventually detect the presence of dark matter 
contained in the inner/intermediate portion of the galaxy. 

To make this comparison easy, the dimensionless quan- 
tities Vmod, o"mod, produccd by the model should be corre- 
lated with the physical parameters. Since the model quan- 
tities are produced in the dimensionless system having 
GMtot = 1 and i?e,tot = 1) they can be related to the 
observables by the use of the simple scaling law: 

y -Ke,tot 

I GMtot 

Cobs — \/ ~5 • O'mod 

y -n-e,tot 

We then produced a set of rotation- and velocity dis- 
persion curves for six type of galaxies having different 
bulge-to-disk mass ratio (77) and concentrations (^), and 
with a given Mtot = 10^^ Mq, and i?e,tot = 10 kpc. In this 
simulation, all the components have an isotropic velocity 
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Fig. 2. Rotation curves and velocity dispersion profiles for Sa galaxies, varying the concentration of the bulge relative to the 
disk components, and the inclination angle i. The two panels below also show the h^, hi parameters profiles as described by 
the model. All the quantities are referred to a galaxy of M — 10^^ Mq, and Rc — lOkpc. The scaling factor is described into 
the text. 



dispersion. In order to limit the number of free parame- 
ters, we stick to the classical case of a r^/'* bulge embedded 
in an exponential disk. The M/ L ratio of the bulge is as- 
sumed to be twice the corresponding value of the disk. 
In the next paragraph we will show the relevance of an 
anisotropic velocity dispersion in modifying the observed 
curves. 

The curves were calculated for the galaxy seen edge- 
on (90°) and at an inclination of 60° and 30°. The choice 
of the rj values has been made trying to reproduce the 
main Hubble stages, from SO to Sd, in agreement with the 
values deduced by Simien & de Vaucouleurs (1986). The 
whole grid of calculation is available on request in elec- 
tronic form. We are able to reproduce the whole sequence 
assuming different slopes or different luminosity laws for 
the bulge and the disk but of course this profusion of data 
cannot be presented here. 

In Fig. H we show the kinematic curves for the case of 
Sa galaxies (here simply defined as galaxies having an = 
0.68). They are arranged with the curves having different 
^ shown in the same panel. 



The relation between the observed curves with that 
presented in the figures is then: 



Vohs — 



Cobs 



/Aftot/[10"MQ] 
i?o,tot/[10kpc] 

/Mtot/[10"MQ~ 

i?e,tot/[10kpc] 



with masses computed in solar mass units and radii in 
Kpc. 

This procedure allows, within some limits, to consider 
the galaxy as a unique structure, with a typical bulge/disk 
ratio and bulge concentration. However, it fits also into 
the idea that it is possible to distinguish dynamically 
the galaxy type. This approach has been faced by Per- 
sic and Salucci (1996) on the basis of a sample of galaxies 
taken from the literature. Working on the superposition 
of nearly 1000 observed gas rotation curves grouped for 
Hubble type, the authors derived a circular rotation curve 
for every given Hubble type. However, since this work was 
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mainly aimed at the detection of dark matter, more care 
is devoted to the outer part of the rotation curves; in ad- 
dition, a poor statistic is available on the early-type disk 
galaxies SO and Sa. 

5.2. The anisotropy in bulges 

The low rotational velocities found in elliptical galax- 
ies with respect to their flattening (Bertola & Capaccioli 
1977, Kormendy & lUingworth 1977) can be interpreted 
in terms of anisotropy in the stellar velocity distribution 
(Binney 1978). The same indications arose from the study 
of the velocity dispersion profiles (Tonry 1983). For sim- 
ilarity with elliptical galaxies, we can expect that some 
anisotropy could be present in bulge dynamics also. On 
the contrary, the presence of anisotropy in the disk com- 
ponent of the galaxies should be less relevant, due to the 
lower value of the velocity dispersion of disks compared 
with bulges or elliptical galaxies. 

The study of velocity anisotropy in disk galaxies may 
give indications on the differences existing between bulges 
and ellipticals, contributing to understand the formation 
process of this two classes of objects. However, due to 
the composite nature of disk galaxies, to measure the 
anisotropy of the velocity distribution from the observ- 
ables (V, a) is much more difficult than in the case of ellip- 
tical galaxies. Working on the same set of data, concerning 
four SO and Sa galaxies, Kormendy & lUingworth (1982) 
deduced the isotropy of their bulges, while Whitmore, Ru- 
bin and Ford (1985) reached the opposite conclusion. Both 
authors used the V/a test, first proposed by Binney (1976) 
for Elliptical galaxies. The fact that the same data may 
be interpreted in such different ways clearly demonstrates 
the ambiguity of this test, when applied to the bulges of 
disk galaxies. 

In a model of galaxy, the shape of the potential de- 
fines, on the galaxy plane, the shape of the rotation- and 
velocity dispersion- profiles. However, the ratio of the cen- 
tral velocity dispersion to the maximum rotational veloc- 
ity is dependent both from the eccentricity e and from 
the anisotropy parameter /3: one may obtain in a model a 
higher V/a ratio increasing the eccentricity and/or de- 
creasing the anisotropy parameter. As a reverse argu- 
ment, we can deduce the galaxy potential by fitting the 
observed rotation curve along the apparent major axis but 
we cannot deduce from this curve the values of e or p. To 
define the values of these parameters, we need more ob- 
servable constraints. 

One of this may arose from the photometry; in partic- 
ular, from the value of the axial ratio obtained by fitting 
the observed isophotes. In such a case, the p parameter 
can be, in principle, derived from the major axis kine- 
matic curves {V and a) only. Such a case is however so 
rarely found in the literature, due to the great errors that 
can derive from the decomposition algorithms, that the 



photometric methods cannot be considered of general ap- 
plication (Seyfert & Scorza 1996). 

A different approach may allow to deduce the ellip- 
ticity and anisotropy with a simultaneous fit of several 
rotation curves taken at different position angles or with 
offsets parallel to the major axis. This kind of data is avail- 
able in the literature for an increasing number of galaxies 
and may be used to derive the anisotropy in their velocity 
distribution. The model presented here can support this 
kind of approach to estimate the amount of anisotropy 
present in an observed galaxy. 

In the following, we shall discuss two points: first, how 
the anisotropy is involved in the stability of a galaxy; sec- 
ond, how much our model is able to distinguish between 
different anisotropies by means of the rotation- and veloc- 
ity dispersion- curves. 




0.2 0.4 O.B O.B 



b/a 



Fig. 3. Axial ratios and anisotropy of our models. Below the 
solid line lies the physically allowed space of parameters, in the 
simple case of a single homogeneous spheroid, and under our 
hypothesis of constant anisotropy. The maximum (3 allowed 
for different values of the axial ratio in the more realistic case 
of two-components galaxies are shown as starred points. The 
labels show the location of the models described in Fig. m 



5.2.1. Anisotropy and stability 

An instability mechanism may be generated by the local 
excess of kinetic energy going in disordered motions with 
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Fig. 4. An example of different kinematic curves along the major axis (left panels) and along an offset parallel to the major 
axis but shifted by 0.5 -Re. The different curves correspond to an SO galaxy with fixed 5, rj values but different {13, e). The label 
of the different curves refers to the pairs of (/3, e) plotted in Fig. 0. 



respect to the local gravitational potential. Its genesis is 
the following: 

Under the hypothesis that the only ordered motion in 
a galaxy be the rotation around the z axis, the z com- 
ponent of the gravitation force can only be equilibred by 
the z component of the velocity dispersion azz, which is 
univocally determined by the distribution of matter and is 
not dependent from the anisotropy parameter f3 . In such 
a galaxy, the radial component of the velocity dispersion 
is higher as the anisotropy increases, enhancing the im- 
portance of the asymmetric drift. As /3 approaches unity, 
the value of aRR tends to infinity, and so the asymmetric 
drift. Following this trend, the asymmetric drift correction 
to the dynamical equation may overcome the gravitational 
term of equation (6), pulling the galaxy outward and de- 
stroying the dynamical equilibrium. As a consequence, 
once fixed in a galaxy the axial ratio of the two compo- 
nents, not every values of /3 are allowed. The "physical" 
range of this parameter is restricted to < /3 < /3max(e), 
with /3max(e) derived by imposing: 

AK'ymdrift < (H) 

in every point of the galaxy. 

Since the circular velocity is not dependent from the 
anisotropy while the correction for asymmetric drift is re- 



lated to it only by means of a scale factor, Eq. can 
be rewritten as 

fie)^^{P^O)<l~P (12) 

so that the dependence from the f3 parameter is easy to 
handle. 

The /(e) can be computed analytically in the limit 
case of homogeneous spheroids only; the result is shown 
in Fig. ^ with a solid line. In the realistic case, with two- 
components and peaked density profiles, /(e) must be 
computed numerically. We estimated the limiting /3max 
from a set of models with bulges of different ellipticities. 
The results are shown in Fig. || (starred symbols), show- 
ing that the homogeneous spheroids curve can be used as 
a starting point for the analysis of the realistic cases. 

Again, we must point out here that this limiting values 
of P are only valid under our hypothesis of (trr == crcfxj), 
and should certainly be recalculated in different kinematic 
conditions. 

5.2.2. Inferring the anisotropy from observations. 

An observer may estimate the anisotropy present in a 
galaxy only if the difference between models with differ- 
ent anisotropy are greater than the observational errors. 
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Table 1. Parameters from the Dynamical Model for NGC 5866. The labels b and d refers to the bulge and disk components. 



scale 


radius 


axial ratio 


Luminosity 


Mass [10"' Mo] 




lb 


fd 


(b/a)i, (b/a)d 


Lt Ijd 


Mb Ma 


i 


35// 


15// 


0.8 0.15 


62 % 38 % 


5.58 1.72 


71° 



In order to give a rough estimate of the diagnostic power 
of the model in determining the anisotropy parameter /3, 
we produced several models in different regions of the di- 
agram of Fig. ^. Looking at the literature, a conservative 
estimate for the observational errors in kinematic curves 
may be established assuming AV = 20 km/sec and Act = 
40 Km/s. 

Curves of a typical model with different anisotropy and 
bulge flattening are shown in Fig. ^. Looking at this figure, 
one can see that the difference in velocities generated by a 
change of anisotropy of 0.1 produces a negligible variation 
along the apparent galaxy major axis. On the contrary, 
significant variations are observable along offset axes, e.g. 
setting the slit parallel to the major axis but shifted by 0.5 
Re, as shown in the same figure. The highest differences 
are visible in the rotation curve (as great as 40 kms~^) 
while the velocity dispersion curves are less sensitive (< 
10 kms^i). 

The sensitivity of the projected kinematic curves to 
variations of P or e is different for galaxies lying in different 
regions of the (/?, e) plane. We can distinguish three typical 
cases: 

1. Almost spherical bulges Under our assumption of 
homogeneous anisotropy, an almost spherical model 
(with bulge having b/a > 0.85 ) does not allow sta- 
ble solutions with significant anisotropy. 

2. Strongly flattened bulges This is the opposite situ- 
ation (b/a < 0.4). In this case, the diagnostic power of 
the model in determining anisotropy is limited. The 
offset curves give only poor additional information 
with respect to that obtained from the major axis data. 
In real galaxies this mean flattening is observed only 
in the late- type spirals, where the disc is preponder- 
ant with respect to the spheroidal component. As a 
consequence, the bulge gives a low contribution to the 
global kinematics. 

3. Intermediate models The higher efficiency of the 
models in distinguishing the anisotropy of a real galaxy 
from observations is found for galaxies with bulges hav- 
ing intermediate fiattening (0.4 < b/a < 0.85). The off- 
set rotation curve is very sensitive to the (3 parameter 
changes, as described above in the example of Fig. ^ 
The presence of an offset spectrum may delete, in this 
case, the uncertainties on the value of (3 derived from 
observations along the major axis. 



6. A test case: NGC 5866 

We tested our model using a real case of a galaxy for which 
a complete set of photometric and kinematic data exists in 
the literature. We chose NGC 5866, for which rotational 
velocity and velocity dispersion profiles are available along 
the major axis and at different offsets from the nucleus 
(Kormendy & lUingworth 1982). This galaxy is an SO seen 
almost edge-on, showing a prominent dust lane on the 
major axis, clearly visible in Fig. ^ and in the photometric 
profiles (Fig. I). 



[kpc 



40 - 



20 



NGC 5868 



-20 - 





































N 



-60 -40 -20 20 40 60 

r [arc see] 



Fig. 5. An image of the galaxy NGC 5866. The approximate 
position and extension of the kinematic profiles (Fisher 1997, 
Kormendy & lUingworth 1982) are shown superimposed on the 
multi-color image of the galaxy (Peletier and Balcell 1997). 



We took the recent kinematic data published by 
Fisher (1997) as well as that published by Kormendy and 
lUingworth (1982). We adopted the photometry made by 
Peletier and Balcell (1997). The full set of data to be re- 
produced by our model include the profiles of velocity dis- 
persion and rotational velocity taken along 7 different axis 
across the galaxy. The position of the axis where kinematic 
data are available are shown in Fig. ^. The and /i4 pa- 
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Fig. 6. Major and minor axis R-band photometric profiles, compared with the two components best-fit model. The dashed lines 
shown the contribution of the bulge and disk to the total surface brightness profiles (solid line) . The crossed circle represent the 
points which have been excluded from the automatic fit due to the effect of the absorption of the dust lane. 



rameters are also published but not for all the positions 
of the sht. 

The effect of the dust is marked in the innermost 5 
arcsec along the major axis. For this reason we neglected 
the regions more heavily obscured by the dust from our 
photometric fit (see Fig. relaxing the fit along the ma- 
jor axis luminosity profile, where the dust absorption is 
more effective. 

The best-fit we adopted arise from a two-component 
(bulge-f disk) model whose parameters are shown in 
Tab. |l]. With these values, we computed by means of 
our model the expected, self-consistent kinematics of this 
galaxy. According to our expectations, the results repro- 
duce most of the observed feature of the kinematics: 

— In general, the trend of the rotation curves and of the 
velocity dispersion profiles agrees everywhere with the 
model, except in the regions obscured by the dust (see 
Fig. ^ . The explanation of this effect may be under- 
stood if the optical thickness is high. In such a case, 
we only observe the stars on the outer edge of the disk, 
whose line-of-sight velocity components are very low. 



— On the other hand, the presence of the dust is only 
slightly affecting the velocity dispersion, as expected in 
a dynamically hot system such as an early-type galaxy. 
Our profiles agree with the observations even in the 
regions obscured by dust. Along the major axis the 
velocity dispersion profile is monotonically decreasing 
from a peak value of 140 kms~^in the central region 
to (T « 74 kms~^at 40". On the major-axis offset posi- 
tions, slightly higher velocity dispersion are observed. 
The minor-axis data indicate that the velocity disper- 
sion rises marginally with increasing z height (from 
« 150 km s^l to « 158 kms^^at 6") and decrease 
at larger radius. The same behavior is present in the 
minor-axis offset spectra. Fisher (1997) proposed that 
this is an observational effect due to the dust obscu- 
ration: our model seem to indicate that this is indeed 
a real effect, due to the dominance of the cold disk 
component at z w 0. 

— Along the major axis, far from the region heavily ob- 
scured by the dust, the model rotation curves flattens 
according to the observed data. Fitting the observed 
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Fig. 7. Rotational velocity profile of stars and gas as compared to the model of the stellar kinematics, for the axis parallel to 
the major axis. The dashed line on the model correspond to the region obscured by the dust lane. In this plot we also show for 
comparison the best results for an anisotropic model with /3 = 0.1 (long-dashed line). 



maximum rotational velocity of 203 kms ^, we find 
a total mass for the galaxy of (7.3 ± 0.5) • 10^° Mq. 
Fairly rapid rotation is seen even at heights of sev- 
eral kpc above the disk plane. At z = 6", r — 40" 
the offset spectrum has a speed of 184 kms"""^. At 
z = 22", equivalent to an height of 2 kpc above the 
disc plane, the rotation curve reaches ~100 kms^^, a 
value comparable with the mean dispersion cr « 150 
km s"""^ observed in the bulge. This is compatible with 
the isotropic model we are proposing, while no model 
with significant amount of anisotropy has been found 
to be able to improve the fit of the data. 



— No rotation is seen along the minor axis (Vmin = 3 ± 3 
km/sec), so we can likely exclude misalignments be- 
tween the different components or a significant bulge 
triaxiality. 

— The fit of the /13 and /14 parameters are marginally 
satisfactory on the major axis, taking into account 
the effects of the dust; but are poorer on the offset 
at 6". However, having considered the deviations from 
the pure Gaussian shape of the line profile greatly im- 
proved the quality of the V,a fit. 
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Fig. 8. Same as Fig. M, but for the velocity dispersion profiles. Note the bad performance of the anisotropic model. 



7. Conclusions 

Our self-consistent, two-components model of disc galaxy 
may reproduce the photometric and the kinematic prop- 
erties of the disc galaxies of different morphological type. 
In particular: 

1. The model is able to produce rotation curves, velocity 
dispersion profiles and ha, curves for disc galaxies 
from SO to Sd, describing the different models in terms 
of relative weight and the concentration of the bulge 
with respect to the disc. A library of model curves, 
projected at different angles of galaxy inclination with 
respect to the sky plane has been generated. The com- 
parison of these models with observations of real galax- 



ies needs only two scaling factors: the total mass and 
the effective radius. 

2. The anisotropy of the velocity field has been studied 
as instability generator inside the galaxy. We showed 
that for each given value of the bulge axial ratio a sta- 
ble stellar system cannot possess a value of anisotropy 
greater than a fixed value. In the extreme case of spher- 
ical systems, only isotropic distributions of velocities 
are allowed, under our assumption of anisotropy ho- 
mogeneous across the galaxy. 

3. The models show a strong diagnostic power in de- 
tecting anisotropy in the bulges of the early-type disc 
galaxies when spectra outside of the apparent major 
axis are available. The best results are obtainable in 
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Fig. 9. Same as in Fig. ^ but for the perpendicular cut shown in Fig. ^ We avoided to display the anisotropic model here. 
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Fig. 10. Same as in Fig. but for the perpendicular cut shown in Fig. ^ 



the case of bulges of intermediate (0.4 < b/a < 0.85) 
flattening, close to the values observed in real galax- 
ies. In this case, an anisotropy change of 15% may give 
variations even of 40-50% in the offset rotation curves, 
as shown in Fig. ^. 
4. The application of the model to edge-on SO galaxy 
NGC 5866 has been presented. It suggests an isotropic 



distribution of velocity and is able to reproduce the 
global galaxy kinematics (V and a profiles). 
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Appendix A. Gauss-Hermite expansion of the line 
shape. 

One should address the question how the deviations from 
the pure Gaussian profile are best quantified. 

Van der Marel & Pranx (1993) represented the line pro- 
file as a sum of orthogonal fimctions in a Gauss-Hermite 
series. The expansion naturally leads to two parameters 
describing deviations from a Gaussian: a parameter /13 
describing asymmetric deviations, and a parameter /14 de- 
scribing symmetric deviations. 

The LOSVDs f{v) can then be described by a Gaus- 
sian plus third- and fourth-order Gauss-Hermite func- 
tions: 

f{v) = Jo exp(-yV2)(l + hsHsiv) + KH^) (Al) 
with y = {v — Vfit)/(Tfit, and where 

Hsiy) = {2V2y^ - 2,V2y)/V& (A2) 

W4(y) = (V-122/' + 3)/%/24 (A3) 

are the standard Gauss-Hermite polynomials, Wj = 
exp{—y'^/2) X 'Hi{y) are the Gauss-Hermite basis func- 
tions, and /i3 and /14 are their amplitudes. Iq is a normal- 
ization constant. 
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